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Abstract 

Void models provide a possible explanation of the "accelerated 
expansion" of the Universe without dark energy. To make the con- 
ventional void models more realistic, we allow the void, an under- 
dense region around us, to be anisotropic and consider an average of 
the distance-redshift relations over the solid angle subtended at the 
observer. We first show that after taking the average of a form of 
the optical scalar equation (distance equation), the effective distance 
equation we obtain coincides with the one for the Lemaitre-Tolman- 
Bondi universe with a Dyer-Roeder-like extension. We then numeri- 
cally solve the equation to compare with observational data of Type 
la supernovae. We find that anisotropy allows smaller size of void and 
larger Q^- 



1 Introduction 

In recent years, reliable distance-redshift relations up to z ~ 1 have been 
revealed by observations of Type la supernovae [H [2] . Comparisons of the 
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data obtained from these observations with the predictions from Friedmann- 
Lemaitre (FL) homogeneous and isotropic universe models suggest that the 
Universe is in a phase of accelerated expansion, and there must exist uniden- 
tified dark energy that causes the acceleration. The nature of dark energy 
has been investigated by many researchers 

On the other hand, the possibilities have been discussed that the distance- 
redshift relations provided by the observations of Type la supernovae may 
be reproduced by taking account of the inhomogeneities in the matter dis- 
tribution and geometry without dark energy. One of such models is the 
so-called void model, in which existence of a large scale underdense region 
(void) is assumed. Actually, an inhomogeneous distribution of matter of 
scale 200 ~ 300/ h Mpc (the Hubble constant is lOOh km/s/Mpc) was sug- 
gested by several authors in galactic surveys, and a void on scales of 
about 250/ h Mpc was found by Blanton et al. [5] in the SDSS commissioning 
data. These results suggest that we might live in a local void with a ra- 
dius of 200 ~ 300//iMpc. Motivated by these observations, Tomita O O [8] 
sought to reproduce the distance-redshift relations using a model in which a 
homogeneous spherical region around us is less dense than the background. 

More recently, Alnes et al. ^ (also, [TOj [11]) showed that there exists a 
model that provides a good fit to the Type la supernova observations (as well 
as other kinds), assuming that the outside of a spherical region is homoge- 
neous and isotropic but the inside of it is described by the Lemaitre-Tolman- 
Bondi (LTB) spherically symmetric dust solution [T^l [131 • Yoo et al. [T5] 
showed that if the arbitrary functions in the LTB solution are determined so 
that the distance-redshift relation in the LTB spacetime coincides with that 
of the ACDM model, the resulting model has a void structure on scales of a 
few Gpc. 

The essential mechanism for these void models to be able to reproduce 
the distance-redshift relations provided by the observations is that the local 
Hubble expansion rate monotonically decreases as the distance increases so 
that the incoming light rays experience effective increase in the expansion 
rate. 

In reality, however, the matter distribution in the Universe is not isotropic, 
and the distance-redshift relation inferred by the observational data should 
be considered as an average over all the directions of sight. This point seems 
to have been largely ignored in the literature. To take into account the ef- 
fect caused by these anisotropics, we in this paper start with the optical 
scalar equation in general inhomogeneous and anisotropic spacetimes, and 
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then take an average of a form of this equation (distance equation) over the 
spheres on the past hght cone. We find that after some estimates, the re- 
suhing equation coincides with the one in the LTB universe with the matter 
density pltb formally replaced by a(z)pLTB, where a{z) is a certain function 
of z. In other words, the equation coincides with the one obtained as a Dyer- 
Roeder-like [TOl [IT] extension of the distance equation in the LTB universe. 
With this effective equation, we numerically investigate light properties of 
void models including effects of anisotropy. 

Mattsson [18] has recently discussed a similar generalization of the Dyer- 
Roeder equation, arguing principally that the original Dyer-Roeder extension 
did not take into account spatial inhomogeneities in expansion rate. To in- 
corporate these inhomogeneities, he replaced the Hubble expansion rate ifpL 
in the Dyer-Roeder equation with Pm{z)Hfi^, with /3m(^) being an arbitrary 
function of z. We do not have to artificially introduce a phenomenologi- 
cal function like l3u{z), since we start with a general inhomogeneous and 
anisotropic setting, which naturally leads to an equation with the effect of 
inhomogeneities in both matter density and expansion rate. Moreover, the 
inhomogeneity and anisotropy we consider is supposed to be made by smooth 
distribution of dark matter, rather than the spiky dumpiness by galaxies as 
assumed in Mattsson [18] and in Dyer and Roeder [161 fTTj . 

This paper is organized as follows. In section [21 we derive our effective 
distance equation by means of averaging. In section [21 we consider a specific 
void model and numerically solve the effective distance equation to search the 
parameter space for the best values that are most likely to explain observa- 
tional data. Section [H is devoted to conclusion and some remarks. Appendix 
[Al provides a review for the original Dyer-Roeder extension. Throughout this 
paper, we use c = 1 as the speed of light. 

2 The effective distance equation for inhomo- 
geneious and anisotropic spacetimes 

In the first subsection below, we transform the optical scalar equation to de- 
rive an equation for the angular-diameter distance Da in an inhomogeneous 
and anisotropic dust universe. We then specialize the equation to the LTB 
spherically symmetric dust universe case for later use. In the second subsec- 
tion, we take a spherical average of the general distance equation to make 
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the equation suitable for comparisons with observational data. 



2.1 Distance equation and its specialization in the LTB 
spacetime 



We start with the optical scalar equation [12], which determines the variations 
of the cross sectional area of a null geodesic congruence. Let us assume that 
the null geodesies in the congruence arrive at the observer through paths 
that are sufficiently far from strong inhomogeneities such as galaxies. This 
allows us to ignore the gravitational lensing effect, and in turn, the shear and 
rotation of the congruence. The optical scalar equation in this case becomes 

^ "^'^ + R,^n^ = ^, (1) 



where A is the cross sectional area of the congruence, the null geodesic 
generator, A the affine parameter, and R^v the Ricci tensor. On the other 
hand, Einstein's equation 

Rt.u - ^Rg^u = SnGT^,, (2) 
where T^^^, is the energy momentum tensor, implies 

R^^k^'k" = SnGT^^k^'k'', (3) 



since k^ is null. With the relation Da oc for the angular diameter 
distance Da, we obtain 

7^^ + 87rGT^.Pr = (4) 
JJa "A 

from Eqs.([T]) and ([3]). The only energy content we consider is dust, so that 
Tfiu = pu^u^ , where p is the energy density for dust, the 4- velocity of 
dust. The frequency of the light observed by a comoving observer at a given 
A with 4-velocity is given by ci;(A) = —u^k^^ (e.g., [SO])- Therefore the 
second term in the left hand side of Eq. (jlj) yields 

T^^Pfc^ = uo{\fp 

= u{Q)\l+zfp, (5) 
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where we have used the definition of redshift z ~ — 1. Thus, Eq.(jlj) 
can be written as 

C}Da + ^-kGpDa = 0, (6) 
where is the second order hnear derivative operator defined by 

1 

We stress that this equation is vahd for any dust universe, as long as the grav- 
itational lensing effect along the geodesies we consider can be ignored. Note 
also that solving this equation requires another equation (geodesic equation) 
to relate A with z, for given p. 

Next, we specialize the last equation to the LTB spherically symmetric 
spacetime case with the observer at the symmetry center. The essential part 
of our task is to convert d/dX to d/dz for the case. 

The LTB metric is given by 

^ _^^2 ^ [^y^'-f dr' + R{t, r fdn\ (8) 
1 + 2L{rj 

where the 'energy function' E{r) is an arbitrary function concerning curva- 
ture. The circumferential radius R{t, r) is the solution for the generalized 
Friedmann equation 

}LR^(^t,r)-^^ = E{r), (9) 

where the 'mass function' M(r) is another arbitrary function, and ' and ' 
denote the derivatives with respect to t and r, respectively. 

Using the analytic function S{x) introduced in Ref.[2T], the solution R 
can be represented simply as 

R{t, r) = {m{r)f'\t - ts{T)fl^S [-2E{t) (^^5^) j , (10) 

where the 'Big-Bang function' ts(r) is a third arbitrary function. The func- 
tion S [x) is the unique solution of the first order ordinary differential equation 

4 , ^ , , dSix) ^0 1 , , 

_(5W+.-iif + 3x-^ = (11) 
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that intersects the x = axis transversely [2T]. This function is non-negative, 
monotonically decreasing, and defined for x < (vr/S)^/^. 

As for the dust matter, the 4-velocity is given by = (1, 0, 0, 0), and the 
energy density by 

M'(r) 

Now, suppose that the hght we observe progresses along a path x^{X) 
for decreasing A with A = corresponding to the observer, so that k'^ = 
—dx^/d\. Then, we obtain 



d dz d 1 duj{\) d 

(13) 



d\ dXdz dX dz 

-1 f dH\ d 



u;(0) \d\^J dz' 

since uj{X) = k^{\) = —dt/d\. On the other hand, the radial geodesic 
equation for t(A) and the null condition, respectively, give us 

dH -R'R fdrV f dtV R'^ fdrV ^ 

= 0. (14) 



dA2 l + 2E\d\J ' \d\J l + 2E\dX 

Using these equations, we immediately have 

d I R' f dtV d 



dX uj{0) R \dX J dz 



= u;{Q)Hj^TB{z){l + zy-^, (15) 

clz 

where we have defined the effective Hubble function Hutb, 

R , , 

^LTB = (16) 

Consequently, Eq.Q yields 

Cl^^Df'Xz) + 4vrGpLTB(^)I^?™^(^) = 0, (17) 
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where 



d 

dz 



d 

dz 



(18) 



and Z)^ represents the angular-diameter distance in the LTB universe. 
Note that this equation is the same as the one for the FL universe (See 
Eg. (162!) in appendix [A]) with the Hubble function replaced by Eq. (fT6!) . (This 
is the reason we call -/^ltb the effective Hubble function.) Note also that the 
effective Hubble function ifLTB(-2) coincides with the so-called longitudinal 
Hubble function H\\ = R'/R'; in contrast, the perpendicular Hubble function 
H_i = R/R plays no role in determining the distance. 



2.2 Averaging of anisotropies 

In this subsection, we take a spherical average of the distance equation 
and compare it with the one for the isotropic spacetime Eq. flTTl) . 

For our purpose, the most convenient is to take a particular coordinate 
system that is based on a foliation of spacetime by past light cones of a 
timelike curve. We employ a spherical null coordinate system (r, z, 6, 0) such 
that T = constant on each light cone, that z = corresponds to the vertex of 
each light cone, that the timelike curve the vertices z = comprise coincides 
with the orbit of matter flow there, and that the angular coordinates {6, 0) 
are constant along each null geodesic that reaches the vertex z = 0. (This 
choice of angular coordinates is possible at least in the domain where the 
null geodesies do not intersect each other.) The radial null coordinate z is 
chosen to give the redshift along each null geodesic labeled by (r, 9, 0). Our 
observation is supposed to be done at the coordinate center 2; = at a certain 
instant r = tq. 

Let Da{z ; 6, 0) be the angular-diameter distance to an object observed in 
the direction {6, 0) with redshift z. This is a function of z, 6, and 0, but we 
will often regard this as a function of z with 6 and being parameters, since 
we mostly consider its variations along the null geodesic, in which situation 
6 and are constant. We will therefore use d/dz instead of d/dz when it 
acts on Da{z ',0,4')- The same rule will be applied to other functions as well. 

Let us define the spherical average of a function f{z ; 6, 0) by 

f{z) ^ {f{z; e, 0)) = ^ 1 /(^; 0, <l>)dn, (19) 
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where dQ = sin 6d6d(f). Then, the spherical average of Eq.([n]) can be written 



as 



{C^Da{z; 6,(1))) + ArtG{p{z; 6, (f))DA{z; 6, 0)) = 0. (20) 
Here, the derivative operator becomes under our coordinates 

2 



CU2(0)(1 + Z)2 \dXdz 



dz d 



H{z-e,m + ^fj-^ 



(21) 



where the general Hubble function H{z ; 6*, 0) is defined by 

"^^-'''^^-zmkw^Tx- '''' 

Note that in the averaged equation fl2Ul) . each term in the left hand side 
depends only on z, as in the LTB equation (fT7|l . We wish to estimate how 
much different the averaged distance function Da{z) = {Da{z ;9,(j))) is from 
an LTB distance function D)^ (z) . To this, we match an LTB solution with 
the given inhomogeneous and anisotropic universe according to the following 
correspondence: 

p(^) = Pltb(^), S{z) = H^tb{z). (23) 

Since the LTB solution possesses two arbitrary physical functions (one of 
the three arbitrary functions corresponds to gauge), it is natural to expect 
that these two conditions are satisfied for a suitable choice of these functions. 
(Strictly speaking, asserting this leads to a variation of the inverse problem 
for the LTB solution, which we consider out of our scope. For a discussion 
of the standard version of the inverse problem, see, e.g., [22] •) 

Keeping this correspondence in mind, we formally rewrite equation (l20i) 
in the form 

{C'Da){z) + AnGaiz)p^T:j,iz)DAiz) = 0, (24) 



where 



_ {p{z;9,<P)Da{z;9,<I))) 
Pltb{z)Da{z} 
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To evaluate the average {C^Da) and the ratio function a{z), let us de- 
compose physical functions as 

p{z- e, 0) = Pltb{z)[1 + S,{z; 6, 0)], (26) 
H{z; e, </.) = i/LTB(^)[l + Sh{z; 9, </.)], (27) 
Da{z; d, (t>) = Da{z)[1 + 5d{z] 9, 0)], (28) 

where (5*) = (* = p,H,OT D). To calculate the quantity C^Da, it is most 
convenient to use the following identity that is confirmed by a straightforward 
computation; 

£^ = (1 + Snm^B + //^tb(1 + znS'^ + ^(4)1^- (29) 
Here, a prime stands for d/dz. Using this formula we can calculate 

C^Da ^C^Da{1 + 5d) 

^{jC'DaKI + 5d) + DaC\1 + 5d) + 2Hl^^{l + Snfil + zYD'aS'i, 
=(1 + 5Hf[{l + 5d)CI^^Da + DaCI^^Sd] 

+ ff^TB(i + + \{^iym + 5d)5:i + daS'u] 

+ 2Hl,^{l + zf{l + 5H?D'^5'^. 

(30) 

(To obtain the second line, Leibniz's rule was applied.) Taking average of 
this equation leads to an expression for (C'^Da), which we, for convenience, 
write in three parts depending on order in S^; 

{C^Da) = CItbDa + 0(2) + 0(3), (31) 

where 0(2) and 0(3) are, respectively, the second and third order term. Brief 
calculations show that 

0{2) = {{5l) + 2{5h6d))CI,^Da 

+ ^ltb(1 + + + {W)D'a (32) 

+ [2{5HCl,^5n)) + //^tb(1 + zf {5'^^ D a, 
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and 



0(3) = {51Sd)CIt.^Da 



Since the perturbation of distance is mainly caused by a perturbation of 
the expansion rate, we may have 6d ~ 5^:/ as a rough estimate. Simphfying 
the second order term according to this approximation and omitting the third 
order term, we have 



(34) 



On the other hand, the function a{z) can be simply calculated as 

«(^) = (\t=^ ( A pLTB{z)[l + 6p{z-d,4>)]DA{z)[l + 6D{z] 9, 

Pltb[Z)-L>a[z) 

= 1 + {5p5d). (35) 

The fact that the emptier the region the light we observe travels, the longer 
the distance from the light source implies that 6p and Sr, have an opposite 
sign each other. Therefore we have 

a{z) < 1 (36) 

with equality being satisfied only in the LTB isotropic universe. 

Furthermore, according to the CMB observation, the perturbation of the 
expansion rate is about \6h\ = \/ {Sjj) ~ 0.1 on scales of 2; ~ 0.1, which 
implies that the second order terms in (C^Da) are much smaller than the 
second order term in 4:nGa{z)puTBDA if \Sp\ ~ 1. We therefore neglect the 
second order terms in (C'^Da) and simply have 

(C'Da) ^ CI^bDa. (37) 

We also note that if ~ | ~ 0.1 and \6p\ ~ 1, then Schwarz's inequality 
implies \{6p6d) \ < \Sp\\6d\ ~ 0.1, or 

aiz) > 0.9. (38) 



10 



To summarize, the distance equation (I24p takes the form 



Clrr^DAiz) + 4nGa{z)p^TB{z)DA{z) = 0, (39) 

where a{z) is a function in the range 

0.9 < a{z) < 1. (40) 

Equation (!39|) is the effective distance equation the averaged distance func- 
tion Da{z) should satisfy. We will use this equation to determine an LTB 
solution and the function a{z), comparing its solution Da{z) with a set of 
observational data. 

Finally, we discuss another possibility where a{z) can be smaller than 0.9. 
Note that if the directions for which the Type la supernovae are observed 
were biased, we would have to restrict the integration region from the whole 
sphere S"^ to part of it, Vi C S^, when taking the spherical average of the 
distance equation fl2Ul) . {Vi is not necessarily a single connected region.) Let 
be the partial average defined by 

{f{z-9,4>))i^ ■ (41) 

Repeating our analysis, the averaged equation 

{C'Da{z; e, 0))i + A7iG{p{z; 6, <j))DA{z; 6, 0))i = 0. (42) 

becomes 

ClDi{z) + AnGai{z)pi{z)Di{z) = 0, (43) 

where Di{z) = {Da)i, Pi{z) = (p)i. The operator is defined with respect 
to the partially averaged Hubble function {H{z;6,(f)))i. Function ai{z) is 
defined in a completely parallel way to the original a{z), which implies 0.9 < 
tti(^) < 1.0. To determine the corresponding LTB spacetime, however, we 
need to consider the same H[z) and p{z) as before with the averaging on the 
whole sphere; H{z) = Hutb{z) and p{z) = putb{z). As a result, we have 

CIt^Mz) + A7TGaiz)p^TB{z)D,{z) = 0, (44) 

where 

a{z)^a,{z)^. (45) 
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(We replaced by Z^ltb a rough estimate.) This is the equation for 
the case of biased distribution of observed hght sources. Apparently, if the 
partial average pi{z) was smaller than the total average p{z) the generalized 
a{z) could be smaller than 0.9. (Actually, any small positive value would 
be possible.) This may be considered as a generalization of the idea of Dyer 
and Roeder [TOl [IT], which is based on a consideration of single light ray, 
whereas the above equation is based on averaging over a partial region in the 
sphere of constant z. Since we do not know whether the directions of SNIa 
are 'biased' or not, in the next section we will consider both possibilities and 
perform our analysis with and without the restriction a{z) > 0.9. 



3 Numerical analysis 

In this section, we numerically seek in a space of control parameters the best 
configuration that can explain observational data of Type la SNe without 
introducing dark energy, and identify the role of anisotropy based on the 
results. 

We employ the same void model considered by Alnes et al. [9J to parametrize 
the LTB solution. That is, we choose the arbitrary functions tsir), M(r) 
and E(r) to be 



tsir) = 0, (46) 

(47) 



M(r) = -H^ 



2 



2Ar 



E{r) = ^HyAM (^1 - tanh , (48) 

where AM is the density contrast, tq the position of the boundary separating 
the inside and the outside of the void, Ar the width of the boundary, and 
Ho the Hubble constant at the center. This void model represents a flat FL 
universe at r ^ ro, since M(r) oc and E{r) ~ there. On the other hand, 
the model represents another FL universe near the center, since M(r) oc 
and E{r) oc there. From the latter fact, we can define the following density 
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parameters at the center: 



= = 1 - fifc. (50) 

These equations allow us to use Vim to parametrize the solution in place of 
AM. We choose the ratio function a{z) as constant, a{z) = a = constant, 
for simplicity. Thus, our parameter space consists of a, Qm, ^o, and Ar. As 
for the observational data to compare, we use the so-called gold data set of 
Riess et al [21], consisting of 157 samples of SNIa. 

To solve the distance equation fl39|) . we need to determine the functions 
-f^LTB(^) and Pltb(^) by solving the null geodesic equations [25] 

dt _ 1 R' dr _ 1 ^/TT2E 

dz~ l + zRi' dz~l + z Ri ' ^ ' 

along the light path {t{z), ^{z)) under initial conditions t(0) = to and r(0) = 
0, with to being the present time. With a solution, the functions Huyb{,z) 
and Pltb(^) can be evaluated according to Pltb(-2) = PuT'B{i{,z)^r{z)) and 
H^r^j,{z) = H^^j,{t{z),r{z)). 

Between the angular-diameter distance Da and the luminosity distance 
Dl holds in general the following relation 



DL{z-e,ct)) = {l + zfDA{z;e,ct)), (52) 

which immediately implies 

DL{z) = {l + zfDA{z). (53) 

This equation allows us to convert the averaged angular-diameter distance 
Da to the averaged luminosity distance -D^ to compare solutions of our dis- 
tance equation with luminosity distance data of SNIa. 

To evaluate the likelihood of given parameters (a, ^2^, "^'o; Ar/ro) to rep- 
resent the right values, we calculate the (s-g-, [I]) 

X^(a, fi™, ro, Ar/ro) = ^tM_ljM)^ ' , (54) 
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where ii{z) is the distance modulus function associated with the numerical 
solution Dl{z), 

fi{z) = 5 log {DLiz)/Mpc) + 25, (55) 
while {zi, fj,°^^{zi), (Tj) is an observational data in the gold data set with Zi be- 



ing the redshift of a supernova, jj, (zi) the corresponding distance modulus, 
CTj the estimated error for fi°^^{zi). The likelihood density is then given by 



P(a, Qrn, ro, Ar/ro) = Ce' 

where C is a normalization constant. 

We searched the parameter region defined by 



-XV2 



(56) 



V = < (a,fim,ro, Ar/ro) 



0.0 < a < 1.0, 0.13 <n„,< 0.40, 
0.15 < ro < 0.42, 0.2 < Ar/ro < 1.0 



(57) 



Numerically solving Eg. (1391) for lattice points in this region, we sought the 
smallest both in whole V and in the region restricted to a > 0.9, V|q,>o.9. 
The best-fit parameters giving the smallest ^ire summarized in Table [TJ 
For comparison, we included the results corresponding to the isotropic void 
model (a = 1). Figures [T] and [2] show confidence regions in selected 2-dim 
parameter spaces, obtained after integrating over the remaining 2-dim space. 
(We determined the confidence regions so that the 100% region corresponds 
to V or V|q,>o.9.) Figure [3] shows residual plots. 



Model 


Isotropic (a = 1) 


Unbiased (0.9 < a) 


Biased (0 < a) 


Ratio parameter: a 


1.0 


0.90 


0.42 


Density parameter: flm 


0.16 


0.17 


0.31 


Transition point: tq 


0.31 


0.26 


0.18 


Transition width: Ar/ro 


0.76 


0.98 


0.47 




175.0 


174.8 


174.1 



Table 1: The best-fit parameters for the isotropic void model (a = 1), the un- 
biased anisotropic model (0.9 < a < 1.0), and the possibly biased anisotropic 
model (0 < a < 1.0). 

Our best-fit value = 174.1 or 174.8 is smaller than Xacdm = ^'^^ 
the ACDM model [21] and Xvoid = 1^6.5 in the void model of Alnes et al. [9]. 
(The value of Alnes et al. is slightly different from our value for a = 1 shown 
in Table [TJ This is perhaps because they sought the best parameters to fit 
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both SNIa and the first peak of the CMB power spectrum, while we only 
consider the SNIa.) 

The most striking feature of a is that smaller value of this parameter 
allows larger and smaller ro (void size). Linear interpolation of the values 
in Table [1] gives us the estimate 

fi^ = -0.26a + 0.42, tq = 0.22a + 0.09. (58) 

The confidence regions shown in Fig{T] also suggest linear relations among 
these parameters consistent with this estimate, although they do not exactly 
coincides with the above estimate. This is because the most likely point in 
the whole parameter space does not necessarily coincide with the most likely 
point in the reduced parameter space that is obtained by integrating the 
likelihood function over uninterested (or unfocused) parameters. As a result, 
the above linear relations do not exactly coincide with the longitudinal lines 
of the confidence contours in FiglH Still, the tendencies are consistent. 

This effect of a may be explained as follows. We know that smaller a 
makes the luminosity distance larger (See Appendix while larger fi^ 
decreases the distance. Smaller void size also decreases the distance, since 
our void is defined as an underdense region, where smaller (possibly nega- 
tive) curvature has an effect of increasing distance. Therefore smaller a can 
compensate larger and smaller void size. Our numerical results seem to 
confirm this explanation. 

The parameter Ar/rg, which is the boundary width of the void (relative 
to its radius), appears to be rather insensitive to the likelihood; the confi- 
dence regions for Ar/ro versus fim, tq, or a shown in Fig{T] do not provide 
reasonably sharp prediction about the choice of Ar/rg. This parameter may 
therefore be considered less important. 

4 Concluding remarks 

Since Type la supernovae are observed in various directions in the sky, the 
distance-redshift relation inferred by a collective data of them should be 
cosidered as an average over the solid angle. On the other hand, the distance- 
redshift relation (distance function) Da{^z) (or Dl{z)) for a given direction 
of sight can be theoretically determined by the distance equation that is 
equivalent to the optical scalar equation. To obtain an effective equation 
for the spherically averaged distance function Da{z) (or Dl{z)), with which 
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the observed distance-redshift relation can directly be compared, we took an 
average of the distance equation, and found that the effective equation we 
want is given by the distance equation for the LTB universe with a Dyer- 
Roeder-like extension. The function a{z) introduced in our effective equation 
represents the degree to which the Universe fluctuates from an isotropic con- 
figuration. Smaller a generally implies larger fluctuation of matter. Our 
numerical computations demonstrated that smaller a allows larger and 
smaller void size. This will provide a useful guideline for further investiga- 
tions of inhomogeneous models. 

In our numerical study, we chose the function a{z) to be constant for 
simplicity, but in reality this function should be a monotonically increasing 
function such that a{z) ^ l{z ^ oo), reflecting the fact that the Universe 
was sufficiently isotropic (and homogeneous) at the last scattering surface. 
In this paper we only considered luminosity distance-redshift data of SNIa 
as obserbational data to compare; for this purpose constant a was suffi- 
cient. However, to include other kind of data like CMB power spectrum it 
is inevitable to have a varying a, resulting in more parameters to search. It 
would still be worth exploring such larger parameter space and compare with 
all available obserbational data including CMB and that of baryon acoustic 
oscillation. This is a future work. 



A The original Dyer-Roeder extension 



Dyer and Roeder [16| llTj considered a simplified model of the Universe in 
which the matter is almost concentrated into objects such as galaxies, and the 
light from a distant object travels through the intergalactic, almost empty 
space. They assume that the geometry is described by the usual Robertson- 
Walker metric 



-dr + a{ty 



+ r\de' + siYi' 



(59) 



where a{t) is the scale factor and K is the curvature index. 

Let (t(A), r(A)) be a null geodesic coming to an observer, with A being the 
affine parameter. Then, the geodesic equation for t(A) and the null condition 
in the FL universe are 



dX^ 



—aa 



dr 
dX 



dx) ^ \dXj 



dr\ 



TV =0' (60) 
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respectively. Using these equations, Eq. (fT^ yields 



d 1 a f dt\^ d 



dX uj{0) a \dX J dz 
a dz 

= a;(0)ifFL(l + ^)'W^, (61) 

where we have used z = ^3^^ — 1 in the second equality, and Hpi^^z) = d/a 
is the Habble parameter for the FL universe. With Eq. fl6T]) . this yields the 
distance equation 

CFi^D^r\z) + 4nGpM^)D^J''\z) = 0, (62) 
where we have defined 



HM^)ii+zr± 



(63) 



TFT ^ 

and Dj^ is the angular-diameter distance in the FL universe and pfl the 
matter density. 

To take into account the effect of lumps. Dyer and Roeder reduced, in 
the above equation, the matter density Pfl(^) on the light path to odrPfl 
by a constant q;dr(0 < «dr < 1)- As a result, the Dyer- Roeder equation is 
given by 

£FLDf''\z) + ArrGaj,KPFL{z)Df''\z) = 0, (64) 

where D^^^\z) is the angular-diameter distance in a lumpy universe with 
parameter odr. In this equation, the value odr = 1 corresponds to the 
homogeneous universe, and the value aDR = to a universe where all the 
matter is concentrated into lumps. 

If the universe is flat {K = 0), we can easily solve Eq. (!64|) under the 
initial conditions ^^^^(0) = and dDf^\o)/dz = I/Hq = l/HpLiO). The 
solution is 

Df'^^z) = , ^ Ul + z)-'' - (1 + z)-'^] , (65) 

^ ^ ^ goV2 5 - 24aDR ' 

. 5 - V25 - 24«DR , 5 + ^25^^24^. 

(61 = -. , &2 = ) . (66) 
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From Eg. 0521) . the corresponding luminosity distance is 



Df'^'iz) = „ f^-^'l [(1 + z)-'^ - (1 + z)-"^] . (67) 
iioV25 - 24aDR 

FigH] shows the relation between the distance modulus Eq. flS^ and the 
redshift z for various odr- As seen from these plots, the distance increases 
when a < 1 as compared to a = 1. (Actually we can show that the distance 
monotonically increases as odr decreases, i.e., dDA{z; aD^) /da^ji < 0, for 
this case.) In other words, the distance takes the maximum at odr = 0, but 
it still insufficient to explain the observation. 
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Figure 1: Confidence regions in the space of two parameters obtained from 
the 4 parameter space V, after integrating over the remaining 2-dim space. 





Figure 2: Confidence regions in the space of two parameters obtained from 
the restricted 4 parameter space V|a>o.9, after integrating over the remaining 
2-dim space. 2I 
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Figure 3: Residual plots of distance modulus for the unbiased best model 
{nm,ro,Ar/ro,a) = (0.17,0.26,0.98,0.90) (red solid line), the biased best 
model (0.31,0.18,0.47,0.42) (orange chain line), and the ACDM model 
{Qrn,^A) = (0.25,0.73) (blue longer-interval dashed line) with the base 
horizontal line being the negative curvature FL model with {Qm, ^k) = 
(0.27,0.73). The purple dotted line, corresponding to (fi^, ro, Ar/ro, a) = 
(0.17,0.26,0.98,1.0), which is the same as the unbiased best model except 
the value of a, is included to show how much a smaller a increases the 
distance. The monotonically decreasing green (shorter-interval) dashed line 
corresponds to the Einstein-de Sitter model {flm = !)• The observational 
data, represented by dots and (error) bars, are taken from the gold samples 
of Riess et al.^. 
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Figure 4: The /i-z relations in the original Dyer-Roeder model with obser- 
vational data. The solid (red), chain (green), and dotted (blue) lines corre- 
spond to odr = 1, C(DR = 0.5, and odr = 0, respectively. The dashed (black) 
line corresponds to the ACDM model. The observational data are the gold 
data [21]. 
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